The Effect of Star Formation on Molecular Clouds in Dwarf 
Irregular Galaxies: IC 10 and NGC 6822 



Glen R. Petitpas and Christine D. Wilson 
Department of Physics and Astronomy 
McMaster University, 1280 Main Street West, Hamilton Ontario, Canada L8S 4M1 

ABSTRACT 

We have observed the 12 CO J =2— 1 and J=3— 2 lines at a few locations in 
the dwarf irregular galaxies IC 10 and NGC 6822 using the James Clerk Maxwell 
Telescope. In addition, we have observed the 13 CO J=2— 1 line for IC 10 and the 
first detection of the 13 CO J=3— 2 transition in a Local Group galaxy. The CO 
line ratios in IC 10 are uniform and are consistent with the average line ratios 
observed in M33 at the la level. These low metallicity environments appear 
to be porous to UV radiation and allow for more efficient heating of molecular 
gas by nearby H II regions. The 12 CO J=3— 2/J=2— 1 ratio for the molecular 
cloud in NGC 6822 is higher than those found for IC 10 and M33 and suggests 
that the 12 CO emission is optically thin in this region. This high line ratio is 
likely the result of its location inside a large H II region with low metallicity and 
low gas content, and requires a hydrogen density greater than 10 4 cm -3 and a 
kinetic temperature greater than 100 K. The 12 CO/ 13 CO J=3— 2 line ratio in 
one of the molecular clouds in IC 10 indicates that the gas must have a rather 
high kinetic temperature of about 100 K. In IC 10 we observe structures on a 
variety of size scales that all appear to be gravitationally bound. This effect 
may help explain the rather high star formation rate in IC 10. 



Subject headings: galaxies: irregular — galaxies: individual (IC 10, NGC 6822) 
— galaxies: ISM — ISM: molecules — Local Group 



1. Introduction 

Dwarf irregular galaxies are rather simple systems, free from the complicated dynamics 
that accompany bars or spiral density waves. This simplicity makes them useful for studying 
the interactions between star formation and the interstellar medium (ISM). Temperature, 
density, and abundance may all affect the star formation process, and thus, determining 
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the physical conditions inside molecular clouds is important for understanding how the 
properties of the molecular gas affect the type and amount of stars formed. Unfortunately, 
extragalactic molecular clouds have low filling factors even within beams as small as 15" and 
so their true temperatures cannot be deduced from the observed peak temperature of the 
spectral line. Instead, observations of many rotational transitions of CO and its isotopomers 
can be combined with radiative transfer models to constrain the physical conditions inside 
these clouds. As it is necessary to observe some of the rarer isotopomers such as 13 CO, 
these techniques have been primarily applied to starburst galaxies {e.g. Tilanus et al. 1991). 

Relatively few dwarf galaxies have been detected in CO and the lines that are observed 
are much weaker than those found in spiral galaxies (Combes 1986; Tacconi & Young 1987; 
Arnault et al. 1988). The weakness of these lines is likely due to a lower CO emissivity 
per unit mass of gas relative to spiral galaxies, which is caused by the low metallicity of 
most irregular galaxies (12 + log(0/H) = 8.16 for IC 10, 8.20 for NGC 6822, Lequeux 
et al. 1979). IC 10 is one of the few dwarf galaxies in which strong CO lines have been 
detected (Becker 1990), with line strengths that are greater than those of the LMC, SMC, 
and NGC 6822 (e.g. Israel et al. 1993; Wilson 1994; Israel 1997). The simplest explanation 
of these strong CO lines is that IC 10 contains large amounts of molecular gas (Wilson & 
Reid 1991). Unfortunately, except for the strong CO lines, there is no additional evidence 
to support this explanation. In fact, IC 10 has a rather low gas-to-dust ratio for a dwarf 
irregular galaxy (Schmidt & Boiler 1993). A more likely answer lies in the star forming 
activity. Massey & Armandroff (1995) found that IC 10 has a global surface density of 
Wolf-Rayet stars three times as high as any other Local Group galaxy, comparable to the 
surface density found within active star forming regions of M33. This high concentration 
of Wolf-Rayet stars reinforces the idea that IC 10 is undergoing a rather strong burst of 
star formation. It has been shown both observationally and theoretically that CO emission 
depends on the conditions in the molecular clouds. Maloney & Black (1988) find that the 
ratio of CO intensity to excitation temperature is constant for ensembles of identical clouds. 
Hence, higher temperatures will result in brighter CO emission. The high star formation 
rates could increase the temperature of the molecular gas and more readily excite the CO 
molecules. Such a process would be able to produce the strong CO lines observed in IC 10. 

In this paper we study the physical properties of molecular clouds in environments 
with different amounts of star forming activity located in two Local Group dwarf irregular 
galaxies, IC 10 and NGC 6822. IC 10 appears optically as a patchy nebulous region 
extending approximately 3' x 4'. It has been classified as a dwarf irregular of the Magellanic 
type (IBm; de Vaucouleurs & Abies 1965). It is rather close to the Galactic plane (b ~ —3°), 
which indicates that there is heavy obscuration at visual wavelengths. It is considered a 
Local Group member based on its distance of 0.82 Mpc (Wilson et al. 1996). The bright 
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H II regions (e.g. Sandage & Tammann 1974) indicate that IC 10 is undergoing massive star 
formation. NGC 6822 is also a dwarf irregular galaxy relatively close to the Galactic plane 
(b ~ —18°). It has a bar-like appearance extending approximately 10' x 20' north-south 
with a prominent "wing" located to the southeast. Its distance of 0.49 Mpc (McGonegal et 
al. 1983) also places it in the Local Group. In §|| we discuss the observations and the data 
reduction process. In §Q we discuss the results found for the individual molecular clouds in 
IC 10 and NGC 6822 based on the line ratios and models. We also compare masses obtained 
from flux measurements with virial masses and discuss the validity of the assumptions made 
to aid in the interpretation of the data. In §^ we discuss the implications of these new data 
and the physical conditions of the molecular clouds for the CO-to-H2 conversion factor. 
The paper is summarized in §[5[ 

2. Observations and Data Reduction 

Four regions centered on giant molecular clouds in IC 10 were observed using the 
James Clerk Maxwell Telescope (JCMT) over the period of 1994 October 20 - 24. The 
regions were selected from both interferometric CO J=l— surveys (Wilson & Reid 1991; 
Wilson 1995) and a single dish 12 CO J=2— 1 and J=l— survey (Becker 1990). All four of 
the regions were observed in 12 CO J=2— 1 with follow-up observations for 2 of these regions 
(containing clouds IC 10:MC 1/2 & IC 10:MC 6/7/8) in 12 CO J=3-2 and 13 CO J =2-1. 
A final 13 CO J=3— 2 observation was obtained at a single position located at the peak of 
the 13 CO J =2-1 emission in IC10:MC 1/2. The half-power beamwidth of the JCMT is 
22" at 230 GHz ( 12 CO J=2-l) and 15" at 345 GHz ( 12 CO J=3-2). Data were obtained 
on a five-point cross offset by 11" for the J=2— 1 data and 8" for the J=3— 2 data. One 
giant molecular cloud in NGC 6822 was observed during service observing at the JCMT on 
the nights of 1996 May 5 ( 12 CO J=3-2) and 1994 June 8 ( 12 CO J=2-l). The cloud was 
observed at a single position in the J =2— 1 transition and a five-point cross offset by 8" in 
the J=3— 2 transition. All observations were obtained using the Dwingeloo Autocorrelation 
Spectrometer. 

The data were reduced using the Bell Labs data reduction package COMB. The data 
were binned to 1.5 km s _1 resolution and first to seventh degree baselines were removed. 
The lines were quite narrow (< 30 km s _1 ) and the spectrometer bandwidth was at least 
450 km s _1 , so the removal of rather high degree baselines for some spectra do not introduce 
significant errors into the final line intensities. Third and seventh degree baseline removal 
on the same data yielded line parameters that agreed to within la. The 12 CO J=3— 2 
data were then convolved to simulate a 22" beam. The spectra for each cloud are shown 
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in Figures |l] and |2|. The data were then binned to the same resolution and shifted to the 
same Vlsr (when necessary). This allowed the calculation of line ratios for each channel 
within the spectral line. The ratios for each channel were averaged to obtain the final line 
ratios. Only those channels with a signal-to-noise ratio of at least two in the line ratio were 
used in the calculation. This method helps reduce errors introduced by the existence of 
poor baselines or weak signals. The uncertainty in the value for the line ratio is taken as 
the standard deviation of the mean. As a final step, the line ratios were scaled to the main 
beam temperature scale using the appropriate values for t^mb (see below). A summary of 
the line ratios is given in Table [I] and a channel by channel histogram of the ratios is shown 
in Figure |. 

The calibration of the IC 10 data was monitored by frequently observing both planets 
and spectral line calibrators. Unfortunately, it was later determined that the beam of 
the telescope was decidedly non-circular during 1994 October (see JCMT Newsletter No. 
4, March 1995), enough to affect our calibration. Although the spectral line calibrators 
showed considerable scatter, with individual measurements differing by as much as ±40% 
from standard spectra, observations of a given source were repeatable to within ~5%. The 
wide scatter of the spectral line calibrators may be understood given the large ellipticity of 
the beam and the range in diameter of the calibrators, from compact carbon stars to more 
extended emission in H II regions. Depending on the size of the source, the measured line 
strength may be larger or smaller with the elliptical beam compared to that measured with 
a circular beam. The main beam efficiencies determined using Saturn were 0.60 ± 0.02 at 
220/230 GHz and 0.36 ± 0.01 at 345 GHz. Both efficiencies are substantially lower than 
normal values. Since Saturn had a diameter of 18" during these observations, it provides 
a reasonable match to the size of the emission regions in IC 10. This assumption was 
confirmed by observing the molecular cloud MC 20 in M33, for which observations obtained 
with the JCMT in a more normal state were available; the 12 CO J=2— 1 line from this 
source was indeed low by about 20%. Thus, we adopt these observed efficiencies to convert 
our data to the main beam temperature scale. For the NGC 6822 data (taken at a different 
dates under normal calibration conditions), observations of CRL 2688 agreed with the 
published values, and so we adopt the normal main beam efficiencies from the JCMT Users 
Guide of 0.72 at 230 GHz and 0.58 at 345 GHz. 

Since the line strengths we are comparing are measured with the same beam diameter, 
using the temperature scale would ensure that our observed line ratios are equal to the 
true radiation temperature ratios. However, conversion to from T£ requires knowledge 
of the forward scattering and spillover (t^fss), which is difficult to measure and was not 
attempted during the observing run. As a result, we do not know how to correct the 
published values of t/fss to take into account the calibration problems discussed above. On 
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the other hand, we do have good values for the main beam efficiencies and so an accurate 
conversion to main beam temperature (T MB ) is possible. Using the T MB scale instead of 
scale only changes the 12 C0 J=3— 2/J=2— 1 line ratio by ~5% under normal calibration 
conditions. For this reason, we will use the main beam temperature scale for the line ratios 
throughout this paper. 

3. Molecular Cloud Properties 

3.1. MC 1/2 and MC 6/7/8 in IC 10 

The clouds MC 1/2 are located 60" (240 pc at a distance of 0.82 Mpc) from the center 
of the brightest H II region in IC 10 (no. Ill, Hodge & Lee 1990; Shostak & Skillman 
1989). This H II region has an Ha luminosity of 4 x 10 38 erg s _1 . Conversely, MC 6/7/8 
is not located near a large H II region; it is 86" (340 pc) away from the edge of even the 
nearest moderate sized one (no. 50, Ha luminosity of ~ 10 37 erg s -1 , Hodge & Lee 1990). 
For these two regions (MC 1/2 and MC 6/7/8), the 12 CO/ 13 CO J=2-l ratios agree very 
well, which suggests similar 12 CO/ 13 CO abundances in these clouds. Although the 12 CO 
J=3— 2/ J=2— 1 ratios also agree within the large uncertainties, the larger ratio of 0.90 is 
found for the cloud nearer the bright H II region, while the lower ratio of 0.66 is found in 
the more quiescent cloud. 

Wilson, Walker, & Thornley (1997) found that the giant molecular cloud NGC 604-2, 
located in the very intense star formation environment of the giant H II region NGC 604 in 
M33, had a 12 CO J=3-2/J=2-l ratio that was higher than average for M33 (1.07 ± 0.17 
compared to 0.64 ± 0.07 for clouds without optical H II regions). They also found that a 
cloud located only 120 pc (31") away from NGC 604 did not show an increase in the 12 CO 
J=3-2/J=2-l ratio. In IC 10, MC 1/2 does exhibit an enhanced 12 CO J=3-2/J=2-l 
line ratio, even though it is further away from the weaker H II region no. 111. This result 
suggests that if a similar line ratio enhancement is taking place as in M33, the ISM in 
the metal-poor dwarf irregular galaxy IC 10 must be more porous to UV radiation. This 
effect could allow the line ratio enhancement to occur at greater distances from weaker H II 
regions than is possible in the higher metallicity spiral galaxy M33 (12 + log(0/H) = 8.48 
in NGC 604, Vilchez et al. 1988). 

The IC 10 data and the NGC 6822 data (see §3.2| ) support the conclusions of Wilson et 
al. (1997) that the 12 CO J=3— 2/J=2— 1 line ratio increases towards the centers of bright 
H II regions. The higher transitions can be more readily excited by higher temperatures or 
higher densities, either of which will produce higher line ratios. It is believed that higher 
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temperatures are present near star-forming regions, but without accurate densities for 
these clouds, we cannot assess the extent to which both of these factors contribute to the 
elevated line ratios. Independent temperature information on these clouds (e.g. from the 
new bolometer array SCUBA on the JCMT) should yield a clearer picture of the processes 
at work. 

For the clouds for which we have observed more than one line ratio, we can use models 
to help constrain the physical conditions of the clouds. We will use the large velocity 
gradient (LVG) model, which assumes that the clouds in the beam have internal velocities 
large enough to Doppler shift the CO emission from one part of the cloud sufficiently to 
prevent reabsorption by other CO molecules in other parts of the cloud (Scoville & Solomon 
1974). We used the LVG code RAD written by Dr. L. Mundy and implemented as part of 
the MIRIAD data reduction package. Model grids were run for kinetic temperatures Tk 
= 10, 15, 20, 30, 50, 100, 200, 300 K, abundance ratios [ 12 CO]/[ 13 CO] = 30, 50, 70, and 
covering a range of densities of log[nH 2 (cm~ 3 )] = 1-6 and 12 CO column densities per unit 
velocity of log[iV( 12 CO)/Aw (cm~ 2 km" 1 s)] = 15 - 20 . To calculate the 12 CO column 
densities from N/Av, we use the full- width half maximum velocity. 

Figure ^ shows one of the two solutions found for the line ratios of IC 10:MC 1/2. 
The three ratios obtained using the JCMT ( 12 CO/ 13 CO J=3-2, 12 CO J=3-2/J=2-l and 
12 CO/ 13 CO J=2— 1) are shown as bands that indicate the la upper and lower limits of the 
line ratio. Solutions are indicated by the locations where the regions of the different ratio 
boundaries overlap. Of the 24 possible models, this only occurred twice for the JCMT line 
ratios of IC 10:MC 1/2, at T K = 100 K and for the abundances [ 12 CO]/[ 13 CO] = 50 and 70. 
Both cases indicate a density of 10 4 — 10 5 cm -3 and column densities per unit velocity of 
about 10 17 ' 7 cm~ 2 (km s -1 ) -1 . The average full-width half maximum of the emission lines 
in IC 10:MC 1/2 is 12.2 ± 1.4 km s _1 , which gives a true column density of 6 x 10 18 cm -2 . 

We have also calculated the 12 CO J=2— 1/J=1— ratio by combining our JCMT 12 CO 
J=2-l data (22" beam) and the IRAM 30-m telescope 12 CO J=l-0 data (21" beam) 
obtained by Becker (1990). The similarity in the beam sizes at these two frequencies allows 
a straightforward combination of the data sets without the use of convolution algorithms. In 
an attempt to constrain the abundance better, we then included the 12 CO J=2— l/J—1— 
line ratio into the LVG model. However, instead of overlapping with one of the other two 
existing solutions, the 12 CO J=2— l/J—1— line ratio provided two additional solutions for 
each abundance at lower density (^(H 2 ) = 10 1,5 — 10 2 ' 5 cm -3 ). There is no region where all 
four line ratios overlap. This result may indicate the existence of separate high and low 
density regions that contribute to the emission. The high density region dominates the 
upper J emission while the lower J emission would originate predominantly from the lower 
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density region (as would be predicted by the Boltzmann equation). We know that Galactic 
molecular clouds are clumpy (e.g. Stutzki & Giisten 1990), and there is every reason to 
believe that extragalactic molecular clouds are clumpy as well. If this is the case we can 
perhaps interpret the two distinct emitting regions as representing the high density clumps 
and the less dense inter-clump regions. Unfortunately, the lack of data with resolution 
better than 7"or 28 pc (Wilson & Reid 1991; Ohta et al. 1992; Wilson 1995) precludes 
further discussion of the degree of substructure within molecular clouds in IC 10. 

With only two line ratios ( 12 CO J=3-2/J=2-l and 12 CO/ 13 CO J=2-l), it is not 
possible to put very strong constraints on the physical conditions in IC 10 MC 6/7/8 
without making some assumptions. For example, if we assume a moderate temperature of 
20 K and a [ 12 CO]/[ 13 CO] abundance of 50 (as found for IC 10:MC 1/2), we can estimate 
upper and lower limits for the densities in the cloud. The results of this model are shown 
in Figure [5[ There is a rather broad range of allowed values for density (10 L6 — 10 4 cm -3 ) 
and for 12 CO column densities per unit velocity (10 165 — 10 18 ' 5 cm -2 (km s -1 ) -1 ). The 
full- width half maximum of the emission lines in MC 6/7/8 is 7.1 ± 1.6 km s -1 , which gives 
true column densities ranging from 2 x 10 17 to 2 x 10 19 cm -2 . If we assume a [ 12 CO]/[ 13 CO] 
abundance of 70 (as is also allowed by the more complete IC 10:MC 1/2 model) the only 
change is that the lower limit on column density increases from 2 x 10 17 cm -2 to 7 x 10 17 
cm -2 , narrowing the allowed range of column density slightly. The range in density and 
the upper limit in column density remain essentially the same. Adopting a higher kinetic 
temperature of 30 K results only in a slightly lower upper limit on H 2 density (10 3,6 cm -3 
instead of 10 4 cm -3 ). Adopting a cooler kinetic temperature of 10 K indicates that H2 
densities greater than 10 4 cm -3 are required to produce the observed 12 CO J=3— 2/J=2— 1 
and 12 CO/ 13 CO J=2-l line ratios. As in IC 10:MC 1/2, the 12 CO J=2-l/J=l-0 line 
ratio for MC 6/7/8 does not overlap with the JCMT line ratios, which again suggests we 
may be seeing a combination of high and low density gas in this region. 



3.2. MC 2 in NGC 6822 

The molecular cloud MC 2 in NGC 6822 is located almost in the very center (within 1") 
of the brightest H II region in NGC 6822. Discovered by Hubble (1925), and designated H V 
by Hodge et al. (1988), this H II region has an Ha luminosity of 4 x 10 38 erg s -1 (Hodge, 
Lee & Kennicutt 1989). The 12 CO J=3-2/J=2-l ratio (1.40 ± 0.07) is significantly 
higher than either of the IC 10 line ratios. The H II region H V is weaker than the giant 
H II region NGC 604 in M33 (Ha luminosity of 5 x 10 39 erg s -1 , Kennicutt 1988) yet we see 
a much higher 12 CO J=3-2/J=2-l line ratio in NGC 6822:MC 2. Again, as in IC 10, we 
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see that in the low metallicity environment of NGC 6822, molecular clouds associated with 
H II regions can attain line ratios higher than in higher metallicity galaxies such as M33. 

High J transitions can be readily excited by high temperatures, high densities, or 
some combination of the two. We would expect the temperature of molecular clouds 
near H II regions to be higher than those unassociated with H II regions. These higher 
temperatures would more readily excite the higher J transitions even in relatively low 
density regions. In the low density, hot gas there will be less CO in the lowest rotational 
levels and hence a lower probability of being reabsorbed {i.e. lower optical depth). Without 
separate constraints on the density, it is difficult to assess the role that density plays in 
producing this high line ratio. However, it can be shown from the virial theorem that for 
clouds with a relatively constant surface pressure the product of density and temperature 
(p T) is approximately constant (e.g. equation (2.8d), McLaughlin & Pudritz 1996). This 
result implies that for virialized clouds, an increase in temperature results in a decrease 
in density. It is likely that this cloud in NGC 6822 is in virial equilibrium (Wilson 1994, 
and next section). It also seems likely in the case of MC 2 that the temperature should be 
higher than in typical interstellar molecular clouds because of its close proximity to the OB 
stars generating the large H II region. We therefore conclude that in MC 2 the high line 
ratio is most likely the result of the high temperatures due to the proximity of the bright 
H II region. 

With only one line ratio, a detailed LVG analysis is not possible. If we assume the 
cloud is in local thermodynamic equilibrium (LTE), we would expect line ratios no greater 
than unity for optically thick gas, while optically thin gas can produce much larger ratios 
(e.g. Walker 1991). The observed line ratio indicates that the 12 CO gas is optically thin in 
MC 2. CO is the most abundant molecule after H 2 and is generally found to be optically 
thick in giant molecular clouds. Optically thin 12 CO is quite rare, and likely the result of 
special physical conditions inside an H II region or possibly in starburst galaxies. 

For gas to be molecular, it must be shielded by dust or have a high enough column 
density to be self-shielding. For molecular gas to collapse into stars, it is believed that the 
column density must exceed some critical value (e.g. Mouschovias & Spitzer 1976; Herbig 
& Terndrup 1986). These values of column density are comparable and are believed to be 
on the order of 10 22 cm -2 for the predominantly molecular clouds in the Galaxy (Solomon 
et al. 1987; Elmegreen 1989; McLaughlin & Pudritz 1996). In starbursts, the gas density 
may vary depending on how long ago the burst of star formation began. Initially, the 
collapse of the molecular clouds could give a higher than average gas density. The density 
should decrease as the gas is used up, new hot stars heat and expand the gas, and stellar 
winds blow gas out of the galaxy. This heating and expansion of the gas in a molecular 
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cloud will decrease the density while leaving the column density unchanged (McLaughlin & 
Pudritz 1996). But, increasing the temperature excites more molecules to higher rotational 
levels, which decreases the opacity in the lower rotational levels. We should then expect 
to see optically thick gas in gravitationally sub-critical molecular clouds unless the clouds 
are heated sufficiently by external sources such as the active star forming regions we see in 
NGC 6822. 

This heating effect may be amplified in the low metallicity environment of NGC 
6822. Lower metallicity suggests a lower dust content, which would reduce the amount 
of shielding by dust from UV radiation from the star forming regions. With this reduced 
shielding, the UV radiation is able to heat the molecular gas more efficiently than is 
possible in higher metallicity (higher dust content) environments, which could contribute to 
the reduced optical depth seen in this galaxy. We might also expect optically thin gas in 
regions where there is simply less molecular gas. There is evidence for a deficit of molecular 
gas in the weak CO line strengths of NGC 6822. We thus attribute the low optical depth 
of the molecular gas in NGC 6822 to the low molecular gas content in this metal-poor 
environment. Our data are insufficient to sort out the extent to which each factor (namely 
low gas content and low metallicity) contribute to the low optical depth. 

Optically thick gas (assuming LTE) can only produce line ratios < 1 whereas optically 
thin gas can produce 12 CO J=3— 2/J=2— 1 line ratios of ^ 2.5 for temperatures greater 
than 100 K (Walker 1991). It is possible that only a fraction of the gas in MC 2 is optically 
thin. For example, the cloud could be a mixture of 10% optically thin gas and 90% optically 
thick gas. Assuming the optically thin CO emission is not blocked by the optically thick 
emission, we would then expect our line ratios to be on the order of 1x0.9 + 2.5x0.1 = 
1.15 for a 100 K cloud. A rigorous calculation of the effect of mixtures of optically thin 
and thick gas is complicated and beyond the scope of this paper. However, this simple 
calculation suggests that MC 2 contains at least ~25% optically thin gas. For simplicity, 
we will assume that the observed high 12 CO J=3— 2/J=2— 1 line ratio indicates that the 
entire cloud is optically thin. 

With only one line ratio for MC 2 ( 12 CO J=3— 2/ J=2— 1), we can only place crude 
limits on the temperature and density. However, we do not need to make any assumptions 
about the abundance ratio, since our data do not include any 13 CO observations. Figure |] 
shows a possible solution for a temperature of 100 K. This temperature seems reasonable, 
considering MC 2 is located within a massive star forming region. If we assume a 
temperature any lower than 50 K, we would require densities greater than 10 6 cm -3 to 
achieve the observed line ratios. While this is not impossible, it seems more likely that this 
line ratio is a result of higher temperatures. For a temperature of 100 K, we can put a lower 
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limit of 10 4 ' 2 cm -3 on density and an upper limit of 10 17 cm" 2 (km s" 1 )" 1 on column density 
per unit velocity. The full-width half maximum velocity of the emission lines in NGC 6822 
is 5.4 ± 0.9 km s" 1 which implies a true upper limit on column density of 5 x 10 17 cm -2 . 

If we increase the temperature above 100 K, we see that lower H2 densities are allowed 
by the models. The observed line ratios can be produced with H 2 densities as low as 100 
cm -3 with temperatures of 200 K, but the CO column densities increase to > 5 x 10 18 
cm" 2 . Assuming a number density ratio of n(CO)/n(H 2 ) of 10 -4 (e.g. Genzel 1992) we 
obtain an H2 column density of about 5 x 10 22 cm -2 . This column density is uncomfortably 
close to the critical column density for molecular cloud collapse of ~ 10 22 cm" 2 . As it seems 
unlikely that we are fortunate enough to observe a cloud that is on the brink of collapse, we 
will use the above value of the CO column density of 5 x 10 17 cm" 2 in further calculations. 



The intensity of the CO emission can be related to the molecular mass of a cloud using 
the equation 



(Wilson & Scoville 1990; Wilson 1995) where S co is the 12 CO J=2-l flux in Jy km s" 1 , R 
is the 12 CO J=2— 1/J—l— line ratio, v is the frequency of the emission (230 GHz for the 
J=2— 1 transition), dyi vc is the distance to the cloud in Mpc, a is the CO-to-H 2 conversion 
factor for the galaxy, and «Gai is the Galactic value (3 ± 1 x 10 20 cm _2 (K km s" 1 )" 1 , Strong 
et al. 1988; Scoville & Sanders 1987). The IC 10 12 CO J=2-l/J=l-0 line ratio was 
calculated by combining Becker's (1990) IRAM 12 CO J=l-0 (21" beam) with our JCMT 
12 CO J=2-l data (22" beam)(see Table 0). For NGC 6822, a 12 CO J=2-l/J=l-0 ratio 
of 0.7 is assumed. The CO-to-H 2 conversion factor (a) is a globally averaged property of 
the galaxy and hence there are uncertainties involved in its use in one specific region of the 
galaxy. Recent studies have suggested that the global CO-to-H 2 conversion factor in IC 10 
may be much higher than in our own Galaxy (a/acai — 40 - 100, Israel 1997; Madden et al. 
1997). However, the conversion factor has been directly measured on smaller scales for the 
regions we are studying by Wilson (1995), and so we adopt her values of a/acai = 2.7 ± 0.5 
for IC 10 and 2.2 ± 0.8 for NGC 6822 and note that this gives the lower limit on the mass 
of the molecular gas. We use 34.6 Jy K" 1 (r? ap = 0.45) and 29.4 Jy K" 1 (r/ ap = 0.53) to 
convert the JCMT data from Kelvins (T£) to Janskys in IC 10 and NGC 6822 respectively 
(Kraus 1986). We assume a coupling efficiency (r] c ) of 0.7 to correct our observed fluxes to 
true fluxes. The virial masses from Becker (1990) were corrected to the new IC 10 distance 
of 0.82 Mpc. The CO-to-H 2 conversion factor is only accurate to within ~ 30% while our 
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Molecular Masses 
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fluxes are typically accurate to about 10%. We therefore adopt a total uncertainty of 40%. 
The results are shown in Table |2|. 

It is interesting that the masses calculated here for the IC 10 regions are larger than 
the virial masses (M v n) found by Becker (1990). Our beam size of 22" covers a diameter of 
88 pc at the distance of IC 10. Although this beam is smaller than the dimensions of the 
clouds in IC 10 as measured by Becker using a comparable (21") beam size (e.g. IC 10:MC 
1/2 measured 107 pc x 98 pc), the IRAM dimensions are not much larger than the beam 
and were not deconvolved from the beam. Thus our single position measurements probably 
detect most of the emission in each region. Given the large uncertainties, the molecular 
mass of all the clouds in IC 10 lie within 2a of the virial masses calculated by Becker (1990). 
If the clouds were not gravitationally bound, the virial masses would be larger than the 
molecular masses. We thus conclude that these quite large structures are gravitationally 
bound. 

The mass of IC 10:MC 1/2 (~ 5 x 10 6 M Q ) indicates that it may be either a large 
giant molecular cloud (GMC, M = 10 5 - 10 6 M , Sanders; Scoville, & Solomon 1985) or 
a small giant molecular association (GMA, M = 10 7 - 10 8 M Q Rand Sz Kulkarni 1990). 
For star formation to occur, we would naturally expect that the molecular clouds in which 
the stars form must be gravitationally bound. All of the GMCs observed in IC 10 appear 
to be gravitationally bound (Wilson & Reid 1991; Wilson 1995). The larger region IC 
10:MC 1/2 also appears to be gravitationally bound. Small GMAs seen in M33 all appear 
unbound by factors of about five except in giant H II regions (Wilson & Scoville 1989; 
1992). In M51, the star formation efficiency is enhanced in the spiral arms where GMAs 
are seen (Rand & Kulkarni 1990). Individual GMCs inside the potential wells of the bound 
GMAs may be colliding at a much higher rate than isolated GMCs, which could increase 
the star formation rate and the star formation efficiency Thus the presence of large bound 
molecular structures in IC 10 may provide a possible explanation for the high star formation 
rate in IC 10. 

For NGC 6822:MC 2, the molecular mass is a factor of 3.4 larger than the virial mass 
calculated by Wilson (1994) using the Owens Valley Millimeter- Wave Interferometer. This 
large mass measurement may be caused by the larger 22" beam overlapping portions of 
the nearby cloud (NGC 6822:MC 1), which is approximately 17" away. Also, the virial 
masses for NGC 6822 are taken from interferometer data, which are insensitive to diffuse 
structures that can be detected by single dishes such as the JCMT. Finally, in determining 
the molecular mass for NGC 6822:MC 2, we assumed a 12 CO J=2-l/J=l-0 value of 0.7. 
Given the rather high 12 CO J=3-2/J=2-l ratio, we might expect this cloud to have an 
unusually high 12 CO J=2— 1/J=1— ratio as well. Increasing the 12 CO J=2— 1/J=1— 
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line ratio to a value of ~1.0 would give molecular and virial mass estimates that agree 
within the uncertainties. 

The CO-to-H2 conversion factor relative to the Galactic value is determined by 
assuming the clouds are virialized and setting the molecular mass M mo \ from equation (JJ) 
equal to the virial mass M vir . This method gives the relation M vir /M mol = a/acai (e.g. 
Wilson 1995). Equation ([]]) assumes that the gas is optically thick and emission comes from 
the surface. However, MC 2 in NGC 6822 appears to be optically thin. For an optically 
thin cloud, the emission from the entire volume of the cloud can escape and we can get a 
much higher flux. This effect would result in an overestimate of the molecular mass found 
using equation ([!]), which in turn would result in an underestimate of a/«Gai- Even if only 
part of the cloud were optically thin, we would still be overestimating molecular masses. 
Because the clouds in NGC 6822 were not fully resolved, the derived virial masses were only 
upper limits and hence the value of a/acai derived by Wilson (1994) is also an upper limit. 
Since NGC 6822:MC 2 is optically thin, we now see that M mo i was also overestimated, 
which would work to balance the previous uncertainty and make the value of a/«Gai found 
by Wilson (1994) more reliable. However, since we do not know the extent of each mass 
overestimate, we cannot make precise corrections to the previously published values of a. 

Given that the 12 CO emission from the cloud NGC 6822:MC 2 is optically thin, we 
have a rare opportunity to calculate the mass of the cloud directly. Using the upper limit on 
column density determined from the LVG analysis, we can estimate the mass by adopting 
a number density ratio of n(CO)/n(H2) = 10~ 4 (Genzel 1992) and multiplying by the area 
of the cloud nab, where a and b are the major and minor axes of the elliptical cloud. It is 
probably most appropriate to use the dimensions of the cloud obtained in the interferometry 
maps (36x < 13 pc, Wilson (1994)) instead of the JCMT beam size. The LVG model 
indicates that the CO column density can be at most 5.3 x 10 17 cm' 2 (at temperatures of 
100 K) to produce the observed line ratio, which leads to a total mass of H2 of < 1.2 x 10 5 
M . This mass is in good agreement with the molecular mass derived using equation ([TJ), 



and agrees with the virial mass to within a factor of two. As discussed in § |3.3| we expect the 
mass derived from single dish flux to be overestimated. Since the virial mass is comparable 
to the mass derived from the column density, it is likely that the molecular cloud MC 2 in 
NGC 6822 is gravitationally bound. 



4. Implications for the CO-to-H 2 Conversion Factor 



It has been proposed that the CO-to-H 2 conversion factor (a) is not only dependent 
on density and brightness temperature (oc ^fp/T h ) but also on the metallicity (Dickman, 
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Snell, & Schloerb 1986; Maloney & Black 1988; Sakamoto 1996). It is possible that the 
lower metallicity in irregular galaxies may change the typical temperatures and densities of 
the molecular clouds they contain. Although this effect may confuse any observed relation 
between a and metallicity, in essence it will simply introduce inter-relationships between p, 
Tb and metallicity. Overall, the expression for a will still be a function of metallicity. 

The line ratios for the clouds in IC 10 agree with those found in M33 by Wilson 
et al. (1997) within the measurement uncertainties (see Table [I]). If we assume that the 
line widths are the same for both transitions, then the line ratios are essentially ratios of 
radiation temperature Tr (e.g. Walker 1991). The radiation temperature is related to the 
cloud excitation temperature and density by 

T R = [J U (T CX ) - J u (T hg )} (1 - e- T ) (2) 

where J U (T) is the Planck function, T bg is the background temperature, and T ex is the 
excitation temperature of the cloud. The optical depth of the cloud is given by 

r D 

t= Kpdr (3) 



o 



where D is the depth of the cloud, p is the density profile of the cloud and K is the 
absorption coefficient. For the optically thin case where r < 1, the 12 CO J=3— 2/J=2— 1 
line ratio becomes 

Tt 2 f 3 ~ 2 iMT^)-MT hg )]T^ 



where / is the beam filling factor of the emission and the superscript 3-2 indicates those 
values pertaining to the J=3— 2 transition. For the optically thick case, r ^> 1, the line 
ratio has the form 

Both these line ratios are a complicated function of the excitation temperature. Thus to 
produce similar line ratios in clouds with different physical conditions would require a 
pathological interdependence of density and temperature to produce the same excitation 
temperature for the low- J CO transitions. We feel that the similarity in the line ratios in 
the clouds in IC 10 and M33 indicates similar physical conditions (namely similar p and T), 
and thus the difference in the measured value of a (Wilson 1995) can be attributed to the 
difference in metallicity of the spiral galaxy M33 and the metal poor dwarf irregulars IC 10 
and NGC 6822. 
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5. Conclusions 

We have observed the 12 CO J=2— 1 and J=3— 2 lines at a few locations in the dwarf 
irregular galaxies IC 10 and NGC 6822. In addition, we have observed the 13 CO J=2— 1 
and J=3— 2 lines for molecular clouds in IC 10. The clouds are located in a variety of 
star formation environments. In general, the clouds in IC 10 show line ratios similar to 
those found in M33 by Wilson et al. (1997). We find that the 12 CO J=3-2/J=2-l ratio 
increases with increasing proximity to H II regions as previously found by Wilson et al. 
(1997) in M33. This result may indicate a general relation that holds true for all galaxies 
(and galaxy types) and is likely the result of a combination of extreme temperature and 
density conditions near star-forming regions. In this low metallicity (and in NGC 6822, low 
CO column density) environment, this effect is more conspicuous than it is in the higher 
metallicity spiral galaxy M33. 

Combining the JCMT line ratios in IC 10:MC 1/2 with large velocity gradient models 
indicates that this region has a kinetic temperature of 100 K, densities of 10 4 — 10 5 cm~ 3 
and a column density of 6 x 10 18 cm -2 . The [ 12 CO]/[ 13 CO] abundance is between 50 and 
70. The inclusion of the IRAM data from Becker (1990) suggests that there are two distinct 
emitting regions, a high density region where the upper J transitions dominate, and a lower 
density region where the lower J transitions dominate. For IC 10:MC 6/7/8, if we assume 
a temperature of 20 K and [ 12 CO]/[ 13 CO] = 50, we obtain densities of 10 16 - 10 4 cm" 3 and 
column densities of 2 x 10 17 — 2 x 10 19 cm~ 2 . 

We find an unusually high 12 CO J=3-2/J=2-l line ratio in NGC 6822 which we 
attribute to optically thin CO emission resulting from the combination of a low molecular 
gas content in a low metallicity environment. An LVG analysis of the 12 CO J=3— 2/J=2— 1 
line ratio in NGC 6822:MC 2 indicates the presence of hot (T > 100 K), optically thin gas. 
The density of this cloud is likely greater than 10 4 ' 2 cm -3 and its column density is less 
than 5 x 10 17 cm -2 . 

The optically thin CO emission in NGC 6822 allows a rare opportunity to calculate the 
mass of molecular gas directly. Mass calculations from two separate techniques determine 
the molecular mass to be 1.2 — 2.1 x 10 5 M Q . This agrees with the virial mass to within 
a factor of ~3. Flux measurements yield masses for the cloud complexes in IC 10 which 
indicate that these regions are gravitationally bound. The masses of the clouds range from 
0.83 — 5.2 x 10 6 M . These regions may represent either large giant molecular clouds 
or small giant molecular associations. The presence of these large gravitationally bound 
structures may help explain the high star formation rate in IC 10. 

The similarity in the line ratios between IC 10 and M33 is most easily understood if 
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the physical conditions in these clouds are also similar, in which case any variation in the 
CO-to-H 2 conversion factor between these two galaxies can be attributed to differences in 
their metallicity. The increase in the 12 CO J=3— 2/J=2— 1 ratio near star forming regions 
may be a result of either higher temperatures, higher densities or some combination of the 
two. Upcoming SCUBA observations will provide temperature measurements which will 
enable us to better determine the physical processes at work. 

The JCMT is operated by the Royal Observatories on behalf of the Particle Physics 
and Astronomy Research Council of the United Kingdom, the Netherlands Organization for 
Scientific Research, and the National Research Council of Canada. This research has been 
supported by a research grant to C. D. W. from NSERC (Canada). 
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Fig. I.— The 12 C0 J =2-1 (solid line), 13 C0 J =1-1 (dotted line), and 12 C0 J=3-2 
(dashed line) spectra observed in 3 regions in IC 10 and NGC 6822. The J=3— 2 data have 
been convolved to match the 22" beam of the J=2— 1 data. The 13 CO J=2— 1 line has been 
scaled up by a factor of three. 

Fig. 2.— The 12 CO J=3-2 (dashed line) and 13 CO J=3-2 (solid line) spectra observed 
towards the region MC 1/2 in IC 10. The spectra were obtained with a 15" beam. The 13 CO 
J=3— 2 spectrum is scaled up by a factor of three. 

Fig. 3. — Channel by channel line ratios (in Tub) for the 3 regions observed. The solid line 
represents the 12 CO J=3— 2/J=2— 1 line ratio (scaled up by a factor of 10), the dashed line 
represents the 12 CO/ 13 CO J=3-2 line ratio, and the dotted line represents the 12 CO/ 13 CO 
J=2— 1 line ratio. Only those channels with a signal to noise ratio greater than two are 
shown. 

Fig. 4. — Solutions for the density and column density for the region IC 10:MC 1/2 using a 
large velocity gradient code. This figure shows the solution for the three JCMT line ratios 
at a kinetic temperature of 100 K and a [ 12 CO]/[ 13 CO] abundance of 50. The solid lines 
indicate the la limits of the 12 CO/ 13 CO J=3-2 line ratio of 5.3 ± 0.4. The dotted lines 
represent the 12 CO J=3— 2/ J=2— 1 line ratio of 0.9 ± 0.11 and the dashed lines represent 
the 12 CO/ 13 CO J=2-l line ratio of 8.2 ± 0.6. Also included is the 12 CO J=2-l/J=l-0 
line ratio of 0.64 ± 0.05 (J— 1—0 data from IRAM, Becker 1990), shown as the dot-dashed 
lines. Solutions are indicated by the shaded regions where the line ratios overlap. In this 
case there are three separate regions where three of the four line ratios overlap which may 
indicate two or more distinct emission regions within the clouds. 

Fig. 5. — Typical LVG solution for the region MC 6/7/8 in IC 10 for a kinetic temperature 
of 20 K and a [ 12 CO]/[ 13 CO] abundance of 50. The dotted lines represent the la limits of the 
12 CO J=3-2/J=2-l line ratio of 0.66 ±0.11 while the dashed lines indicate the 12 CO/ 13 CO 
J=2— 1 line ratio of 8.3 ± 1.4. The shaded region indicates the solutions for the JCMT data. 
Also included is the 12 CO J=2-l/J=l-0 line ratio of 0.50 ±0.06 ( 12 CO J=l-0 data from 
Becker 1990). The 12 CO J=2-l/J=l-0 line ratio does not overlap with either of the JCMT 
line ratios. 

Fig. 6. — This figure shows one possible temperature (100 K) solution for the region MC 2 
in NGC 6822. The dotted lines show the la limits of the 12 CO J=3-2/J=2-l line ratio 
of 1.4 ± 0.07. With only a single line ratio, it is only possible to place lower limits on the 
temperature (> 50 K) and density (> 10 4 cm -3 ). 
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Table 1. CO Line Ratios of Molecular Clouds in IC 10 and NGC 6822 



Cloud a 


a(1950) 


5(1950) 


12 CO J=2-l 


12 CO J=3-2 12 CO J=3-2 


12 CO J=2-l 


13CO J=2-l 


12CO J=2-l 13 CO J=3-2 


12CO J=l-0 


IC 10:MC 1/2 


00:17:46.2 


+59:00:28 


8.2 ± 0.6 


0.90 ± 0.11 5.3 ± 0.4 


0.64 ± 0.05 


IC 10:R2 


00:17:34.7 


+59:02:05 






1.18 ± 0.1 


IC 10:MC 4 


00:17:37.5 


+59:04:08 






0.65 ± 0.06 


IC 10:MC 6/7/8 


00:17:40.1 


+59:04:38 


8.3 ± 1.4 


0.66 + 0.11 


0.50 ± 0.06 


NGC 6822:MC 2 


19:42:03.3 


-14:50:27 




1.54 ± 0.08 




M33 average 13 






7.2 ± 1.1 


0.69 ± 0.15 


0.67 ± 0.19 



All data obtained at JCMT 15 m, except for the 12 CO J=l— data which are taken from Becker 
(1990). All line ratios refer to a 22" beam, except for the 12 CO/ 13 CO J=3— 2 ratio which is for a 
15" beam. The line ratios are the average of all channels with a signal to noise greater than two 
(see Figure |3|) and the uncertainties are the standard deviation of the mean. 

a Cloud designations from Wilson & Reid 1991, Wilson 1994, Wilson 1995 except R2 (Becker 
1990). 

b From Wilson et al. 1997, Thornley & Wilson 1994. 
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Table 2. Masses of Molecular Clouds in IC 10 and NGC 6822 



Cloud 


<S'cO(J=2-l) 

(K km s^XTX) 


Molecular Mass a 
(M©) 


Virial Mass b 
(M©) 


IC 10:MC 1/2 


9.16 


5.2 x 10 6 


2.5 x 10 6 


IC 10:R2 


2.67 


0.83 x 10 6 


0.35 x 10 6 


IC 10:MC 4 


2.37 


1.3 x 10 6 


0.54 x 10 6 


IC 10:MC 6/7/8 


2.35 


1.7 x 10 6 


0.66 x 10 6 


NGC 6822:MC 2 C 


1.66 


0.21 x 10 6 


< 6.3 x 10 4 



a A CO-to-H2 conversion factor (a/acai) of 2.7 and 2.2 has been 
assumed for IC 10 and NGC 6822 respectively. 

b Virial masses taken from Becker (1990) except NGC 6822:MC 2 
(Wilson 1994). 

C A 12 CO J=2-l/J=l-0 ratio of 0.7 has been assumed. 
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